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Abstract 

In the paper, we apply the kx factorization approach to deal with the B — ► K transition form 
factor F+~^ K (q 2 ) in the large recoil regions. The B-meson wave functions ^b and ^>b that include 
the three-particle Fock states' contributions are adopted to give a consistent PQCD analysis of 
the form factor up to 0(l/m 2 ). It has been found that both the wave functions VP b and ^>b can 
give sizable contributions to the form factor and should be kept for a better understanding of the 
B meson decays. Then the contributions from different twist structures of the kaon wavefunction 
are discussed, including the 5C/^(3)-breaking effects. A sizable contribution from the twist-3 wave 
function is found, whose model dependence is discussed by taking two group of parameters 
that are determined by different distribution amplitude moments obtained in the literature. It is 
also shown that F^ K (0) = 0.30 ± 0.04 and [F^ K (0)/ F^ n (0)} = 1.13 ± 0.02, which are more 
reasonable and consistent with the light-cone sum rule results in the large recoil regions. 
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I. INTRODUCTION 



A study on the heavy-to-light exclusive processes plays a complementary role in deter- 
mination of the fundamental parameters of the standard model and in developing the QCD 
theory. And there is an increasing demand for more reliable QCD calculations of the heavy- 
to-light form factors. We have done a consistent analysis of the B — > tx transition form factor 
in Refs. l|, |2J, which shows that the results from the PQCD approach, the lattice QCD ap- 
proach and the QCD light cone sum rules (LCSRs) are complementary to each other and 
by combining the results of those three approaches, one can obtain an understanding of the 
B — > tx transition form factor in the whole physical regions. It is argued that by applying the 
kr factorization approach , where the transverse momentum dependence for both 

the hard scattering part and the non-perturbative wavefunction, the Sudakov effects and 
the threshold effects are included, one can regulate the endpoint singularity from the hard 
scattering part effectively and derive a more reliable PQCD result of the B-meson decays. 
Furthermore, by applying the B-meson wave functions up to next-to-leading Fock state, 
Ref.jsJ calculated the B — > tx transition form factor up to 0(1/ ml) and also made some 
discussions on the reasonable regions for the two phenomenological B-meson wave function 
parameters A and 8, where A is the effective mass of B meson that determines the B-meson's 
leading Fock state behavior and 8 is a typical parameter that determines the broadness of 
the B-meson transverse distribution. Since both pion and kaon are pseudo-scalar mesons, it 
will be interesting to give a consistent PQCD analysis of the B — > K transition form factor 
up to order 0(1/ ml) based on the results of B — > tx transition form factor. 

re B — > K transition form factor has been studied under several 
Ki 11]. A PQCD calculation has been done in Ref.[?], which can be 



In the literature, t 

~k 0, y, 



approaches 



roughly treated as a leading-order estimation 0(1/ m^) since some of the power suppressed 
terms both in the hard scattering amplitude and the B-meson wave function have been 
neglected. The B —> K transition form factor has also been analyzed by several groups 



under the QCD LCSR approach 



ll| ], where some extra treatments to the correlation 



function either from the B-meson side or from the kaonic side are adopted to improve their 
LCSR estimations. New sum rule for the B — > K form factor is derived by expanding the 
correlation function near the light cone in terms of B-meson distributions [8(, in which the 
contributions of the quark-antiquark and quark- ant iquark-gluon components in the B-meson 
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are taken into account. While in Ref.Q] an improved LCSR approach that had been raised in 
Ref. [1^1 was adopted to eliminate the contributions from the most uncertain kaonic twist-3 
wave functions and to enhance the reliability of sum rule calculations of the B* — > K form 

n 

factor. A systematic QCD LCSR calculation has been done in Ref.jlO] by including the one- 
loop radiative corrections to the kaonic twist-2 and twist-3 contributions, and the kaonic 
leading-order twist-4 corrections. It can be found that the main uncertainties in estimation 
of the B —>■ K transition form factor come from the B-meson wave function and the kaonic 
twist-2 and twist-3 wave functions. 

In doing the PQCD calculations on the B-meson decays, an important issue is whether 
we need to take both the two wavefunctions $ b arid into consideration or simply \1> b 
is enough? By taking the frequently used first type definition of \I/ R = ^ B t^ B and ^b = 
fli*. where n are dehned in Ref.fl, it can be found that QQ their distribute 
amplitudes have quite different endpoint behaviors even under the Wandzura-Wilczek (WW) 
approximation 16[ , such difference may be strongly enhanced by the hard scattering kernel. 
For example, the ratio between the contributions of and ^ b is about (—70%) 
for the B — > 7r from factor in the large recoil regions. So the contribution from \1/ b under 
the above definition can not be neglected and it is needed to suppress the big contribution 
coming from $b so as to obtain a reasonable total contributions. To derive more accurate 
estimation, Ref. [tJ raised the second type definition of ^> b = ^b and ^ r = (^b~ ^b)> under 
which the contribution from is of order 0(l/m&) to that of ^ b 18| . For convenience, in 
the following, we shall adopt the second type definition of ^ b and ^/b to do our calculation. 
Then one may ask is it enough to give a 0(1/ ml) estimation with ^b and ^b under 
the WW approximation ? As has been pointed out in Ref.[19j, the 3-particle Fock states' 
contributions to the B-meson wave function can be estimated by attaching an extra gluon to 
the internal off-shell quark line, and then (1/m^) power suppression is induced in comparison 
to that of the WW-part B-meson wave functions. Recently, the B-meson light-cone wave 
functions have been investigated up to next-to-leading order in Fock state expansion in the 
heavy quark limit It was shown that by using the relations between 2- and 3- particle 
wave functions derived from the QCD equations of motion and the heavy-quark symmetry, 
one can give a constraint on the transverse momentum dependence of the B-meson wave 
function, whose distribution tends to be a hyperbola-like curve other than a simple delta 
function that is derived under the WW approximation. These results provide us a chance 
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to give a consistent PQCD analysis of the B — > K form factor up to order 0(1/ ml). 

Another issue we need to be more careful is about the kaonic wave functions. The 
distribution amplitude (DA) for the twist-2 wave function f y has been deeply studied, e.g. 



by the light-front quark model 
calculation 



201 ] . the LCSR approach 



21 



22 



23, 



24j and the lattice 



25l . |26I ] and etc. In Ref.[2lj], the QCD sum rule for the diagonal correlation 



function of local and nonlocal axial-vector currents is used, in which the contributions of 
condensates up to dimension six and the (9(a s )-corrections to the quark-condensate term 
are taken into account. The first Gegenbauer moment af(lGeV) of the twist-2 DA derived 
there, i.e. af(lGeV) = 0.05±0.02, is consistent with that of the lattice calculations 25l. l26|. 
so we shall constrain of (lGeV) within this range when constructing a model for \1/ k- As 
for the twist-3 wave function the calculations of it has more uncertainty than that for 
the leading twist, e.g. its DA moments in Refs. 24], |27|, |28| are quite different from each 
other, where the DA moments in Refs. [24]. I27] are derived by using the QCD light-cone sum 
rules and the moments in Ref. 28|] are derived based on the effective chiral action from the 
instanton. Under the PQCD approach, according to our experience on the B — > 7r transition 
form factor [lj and the pion electro- magnetic form factor 29[, it can be found that for a 
twist-3 wave function with a better endpoint behavior other than the asymptotic one, the 
twist-3's contributions are indeed power suppressed to the leading twist's contribution that 
favor the conventional power counting rules. In the present paper, we shall adopt two groups 



of DA moments 



24 



281 ] together with the Brodsky- Huang-Lepage (BHL) prescription 30] 



to construct a model for ty p , and then make a discussion on its uncertainty to the B — > K 
transition form factor. The S'f//(3)-breaking effects shall also be included for constructing 
the kaonic wave functions. 

The purpose of the paper is to reexamine the B —>■ K transition form factor in the 
PQCD k T factorization approach up to 0(1/ ml). Under the k T factorization approach, 
the full transverse momentum dependence (/^-dependence) for both the hard scattering 
part and the non-perturbative wave function, the Sudakov effects and the threshold effects 
are included to cure the endpoint singularity. Furthermore, we shall analyze the power 
suppressed contributions from both the wave functions and the hard scattering amplitude 
and then give a consistent analysis of the form factor up to 0(1/ ml), which have not been 
considered in the literature. In section II, we give the calculated technology for the form 
factor in the large recoil regions. Also we present the model wave functions of the kaon 
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with better endpoint behavior in the same section, which are constructed based on the BHL 



prescription [30] and the DA moments obtained in Ref. 



24 



28j |. In section III, we give our 



numerical results. Conclusion and a brief summary are presented in the final section. 



II. CALCULATION TECHNOLOGY FOR THE B -> K TRANSITION FORM 
FACTOR 



The B — ► K transition form factors F^ K (q 2 ) and F^^ K (q 2 ) are defined as follows: 
(K(P K )\s lfl b\B(P B )) 



Ml- Ml 



(Pb + Pk),- - 2 K % 
q z 



M 2 — M 2 
F^ K {q 2 ) + MB q2 MK q,Fr K {q 2 ) 

(1) 

where F^~' K (0) should be equal to F^~^ K (0) so as to cancel the poles at q 2 = 0. The 
amplitude for the B — > K transition form factor can be factorized into the convolution 
of the wave functions for the respective hadrons with the hard-scattering amplitude. In 
the large recoil regions, the B — > K transition form factor is dominated by a single gluon 
exchange in the lowest order. In Ref.jlj], we have done a consistent analysis of the B — > 7r 
transition form factor within the kr factorization approach, where the power suppressed 
terms up to 0{l/m 2 ) have been kept explicitly in the hard scattering amplitude. The 
interesting reader may refer to Ref.{l] for more details 1 . More specifically, for the present 
case, one needs to know the momentum projection for the matrix element of the kaon and 
B meson in deriving the hard scattering amplitude. By keeping the transverse momentum 
dependence in the wave function, the momentum projection for the matrix element of the 
kaon has the following form, 



'5 = XK 7b k ^)-^s ki) - «v I n"n v — - p" aK Q ±J — 



M5, = J r\ lis *k(x, kj.)-^ 75 ( %(x, k ± ) - i<v [ n»n v °\ " 7 -f 

(2) 

where fx is the kaon decay constant and fix is the phenominological parameter /i^ = 
M K /(m s + m u ), which is a scale characterized by the chiral perturbation theory. ^k{x, kj_) 
is the twist-2 wave function, ^ p (x, kj_) and kj_) are twist-3 wave functions, respectively. 
^' a (x, kj_) = d^ a (x, k±)/dx, n = (y/2, 0, 0j_) and n = (0, a/2, 0j_) are two null vectors that 



Three typo errors are found in Ref. 1J, i.e. in Eq.(3) £^S£b should be changed to £m+Mm_ > in Eq ^ thc 
factor [3 — r\ — xif\ should be changed to [3 — r\ + xrf\ and in Eq.(7) y should be changed to rj. 
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point to the plus and the minus directions, respectively. And the momentum projection for 



the matrix element of the B meson can be written as 
if B \j) B + M B 



17 



3l|: 







L) + ^(e, ijl) - A(e, w w 



75 



(3) 



a/3 



where £ = is the momentum fraction for the light spectator quark in the B meson 
and A(f,lj_) = M B fH , LJ - 1 ± )). The four-component ^ in Eq.® is 

defined through, = Z M — @lz^±i nM ) with / = (^=, lj_). By including the Sudakov 
form factors and the threshold resummation effects, one can obtain the formulae for the 
B — > K transition form factors F^~* K (q 2 ) and FP~^ K (q 2 ) in the transverse configuration 
6-space, which can be simply obtained from Ref.flJ by changing the pion wave functions 
to the present case of kaon and by changing \1/ B and to the second type definition as 
described in the INTRODUCTION. 

In PQCD approach, the parton transverse momenta are not negligible around the 
endpoint region. The relevant Sudakov factors from both k± and the threshold resummation 



321 ] can cure the endpoint singularity which makes the calculation of the hard amplitudes 



infrared safe, and then the main contribution comes from the perturbative region. Also it 
is necessary to keep the transverse momentum dependence in the wave functions to derive 
a more reliable estimation in PQCD. In principle, the Bethe-Salpeter formalism [33!] and 



the discretized light cone quantization approach 34j could determine the hadronic wave 
functions, but in practice there are many difficulties in getting the exact wave functions at 
present. The BHL prescription [3^], which connects the equal-time wave function in the rest 
frame and the wave function in the infinite momentum frame, provides a useful way to use 
the approximate bound state solution of a hadron in terms of the quark model as the starting 
point for modeling the hadronic wave function. So in the present paper, we will adopt the 
BHL prescription for constructing the kaonic wave functions. While for the 5-meson wave 
function, they have been investigated up to next-to-leading order in Fock state expansion in 



the heavy quark limit in Ref. 2J, which sha 
A simple model has been raised in Ref. 



1 be adopted to do our discussions. 
2J for the B- meson wave functions ^ B and 
which keep the main features caused by the 3-particle Fock states and whose transverse 
momentum dependence are still the like-function of the off-shell energy of the valence quarks 
but shall broaden the transverse momentum dependence under the WW approximation to 
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a certain degree. And in the compact parameter 6 B -space, it reads Q 



and 



167r ) — — exp 



(167T 3 



exp 



^0 , 



T[6\J S -x[k] + (1 - S)T[2 - 5) J!_,[«] - 



l-<5 



(r[5] j,_ x [k] + (i - 5)r[2 - <5] Jx_i[K]) (| 



1-5 



(4) 



(5) 



where Wo = 2A/3, £ = A/M B and « = #(2£ - £)\/£(2£ - £)M B b B . The factor (16vr d ) is 
introduced to ensure their Fourier transformation, i.e. ^b(£, kj_), satisfy the normalization, 
/ 3 X ^g(£, kj_) = 1. It can be found that both and have the same transverse 
momentum dependence and only two phenomenological parameters A and 5 are introduced. 
A is the effective mass of B meson, A = M B — m^, which determines the B meson's leading 
Fock state behavior. 5 is a typical parameter that determines the broadness of the B- 
meson transverse distribution in comparison to the WW-like one. The WW-like B-meson 
wave functions in the compact parameter & B -space can be found in Ref. And a direct 
comparison shows that when 8 — > 1, the transverse momentum dependence of the B-meson 
wave function in Eqs.( H|5|) returns to a simple ^-function, which is the same as that of the 
B-meson wave function under the Wandzura-Wilczek approximation [^j]. According 
to the definitions, we have y B (£,b B ) = Vt(£,b B ), ^(f , &b) = V%(£,b B ) - b B ) and 

A(£,6 B ) = -M B J* d£'* B (£',& B ). 

Next, we construct the kaonic twist-2 wave function based on its first Gegenbauer moment 



af and on the BHL prescription 



the 



25, 



ight-front quark model 



30| • The first Gegenbauer moment af has been studied by 



20] , the LCSR approach 2jJ, |22j, |24J and the lattice calculation 



261 ] and etc. In Ref.[2jJ, the QCD sum rule for the diagonal correlation function of 
local and nonlocal axial-vector currents is used, in which the contributions of condensates 
up to dimension six and the (9(a s )-corrections to the quark-condensate term are taken into 



account. The moments derived there are close to that of the lattice calculation [25|, [26], so 
we shall take af(lGeV) = 0.05 ± 0.02 to determine the model wave function ^k- Based on 
the BHL prescription, we take the twist-2 wave function of kaon as 



*A-(x,kj 



[l + B K C 3 1 /2 {2x- 1)] x 



A 



K 



exp 



K 



kl + ml kl + ml" 



-3 2 ' " " 



x 



1 — x 



(6) 



x(l — x 

where q = u, d, C f (1 — 2x) is the Gegenbauer polynomial. In comparison to the pion wave 
function (e.g. [36|), it can be found that the SUf(3) symmetry is broken by a non-zero B K 



and by the mass difference between the s quark and u (or d) quark in the exponential factor. 
The SUf(3) symmetry breaking in the lepton decays of heavy pseudoscalar mesons and in 



the semileptonic decays of mesons have been studied in Ref.[37j. For definiteness, we take 
the conventional values for the constitute quark masses: m q = 0.30GeV and m s = 0.45GeV. 
The parameters A%, Bk and (3 k can be determined by the value of of together with the 
normalization condition: 

[ l dxf ^v^( x , k± ) = l (7) 

JO Jkl<)J% 167T 3 

and the constraint (k^ 2 « (k-lVJ 2 = 0.350GeV [38], where the average value of the 
transverse momentum square is defined as 

The parameter /io in the model wave function stands for some hadronic scale that is of order 
0(1 GeV). For clarity, we set /i = 1 GeV. The DA <p K (x,fi ) is defined as (f)fc(x, Un) 



J dxd 2 ~k±_ 


kil 


\*k{x, k ± )| 2 


J dxd 2 hj_ 


* K (x,k ± )\ 2 



h 2 ± <nl kj_). The first Gegenbauer moment af-(/i ) of Refs. 21, 22, 24] can be 

defined as 

Kt \ = Io dx<p K {l - x, /i )C 1 3/2 (2x - 1) 
ai tidx6x(l - x)[Cl /2 (2x - l)] 2 ' 

where 0x(l — x.^lq) other than 4>k(x, /io) should be adopted, since in Refs.[2ll. |22j, |24j x 
stands for the momentum fraction of s-quark in the kaon (K), while in the present paper 
we take x as the momentum fraction of the light g-(anti) quark in the kaon (K) 2 . Based on 
the above discussions, we can obtain the values for Ak, Bk and (3 k- 

A K = 2.71 x K^GeV" 1 , B K = [0.116 - 0.9af (#,)] , (3 K = 0.877GeV _1 , (9) 

where the values of Ak and (3k are almost constant, i.e. their changes (5Ak/Ak) and 
(SPk/Pk) are less than 0.001 by varying ctf(/i ) within the range of [0.03,0.07]. More 
specifically, for the case of af (/i ) = 0.05, we have 

A K = 2.71 x 10 2 GeV _1 , B K = 0.071 , (3 K = 0.877GeV _1 . 

As will be seen that the contributions from twist-3 wave function ^ a (x,k±) is less im- 
portant in comparison to that of ^k{x, k±) and ^ p (x, k±), which is similar to the case of 



2 In the literature, there are some ambiguities in use of <pK,p(x, /io) or (f>K,p(l — /to) m connection to the 
hard scattering part. This will cause errors when the iS'C//(3)-symmetry is broken. 
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FIG. 1: Kaon (j) p {x, /io) with its parameters determined by the two groups of DA moments 



24, 



The solid line and the dashed line are for (j>h(x, /xq) and 4> 2 (x, //q) respectively. For comparison, the 



2J and 



|] respectively. The 



big dotted line and the dash-dot line are for <j) s p r {x, no) 
dotted line is the asymptotic behavior of </>p S (x, oo) = 1. 

B — > 7r transition from factor So basing on the BHL prescription, we directly take the 
twist-3 wave function of kaon as 



A CT exp 



-0 



K 



'k\+m 2 q 

x 



+ 



k\ + m 2 s 
1 — x 



(10) 



where A a can be determined by its normalization condition, i.e. A a = 1.36 x 10 3 GeV _1 . 

As for the twist-3 wave function ^ p (x, k±), its DA's asymptotic behavior, 4> p s (x, oo) = 1, 
so its endpoint singularity is much more serious. Then the transverse momentum dependence 
of ty p (x, kj_) is much more important than that of ^k(x, k±) and ty<j(x, k±) in order to cure 
the endpoint singularity. One can construct ^ p (x, k±) in the following form, 



$? p (x, kj 



[1 + B p C\ /2 {2x - 1) + C p C L 2 '\2x - 1)] x 



1/2, 



A~, 



x(l 



x) 



exp 



-0: 



'k\ + m 2 k 2 , + ml 



K 



x 



+ 



1 — X 



(11) 



where x stands for the light quark g's momentum fraction, C^ 2 {2x — 1) and C\ (2x — 1) 
are Gegenbauer polynomials and the coefficients A p , B p and C p can be determined by its 
DA moments. The DA (p p (x,/io) is denned as (j) p (x,/io) = f k 2 <fJ 2 ^^^(x, kj_). 

To discuss the uncertainty caused by ty p , we take two groups of DA moments that have 



been obtained in Refs. 



24 



281 ] to determine the coefficients A p , B p and C p , where the mo- 



Ref. 



ments in Ref. 24j are derived by using the QCD light-cone sum rules and the moments in 



28j are derived based on the effective chiral action from the instanton: 



Group 1 [23] : (x°) 



0\K 
P 



Group 2 [27] : (x°)f 



x 1 )* = 0.06124, (x 2 )J = 0.36757 



(x 1 )* = 0.00678, (x 2 )« = 0.35162. 



2\K 
P 



(12) 
(13) 



Here the moments are defined as (x 1 )^ = J dx(2x — 



24 



)V P (i - x,fio) with i = (0,1,2). It 
281 ] are for P (1 — x, /i ) other than 



should be noted that the moments defined in Ref. 
<fip(x,no), since in these references x stands for the momentum fraction of s-quark in the 
kaon (K), while in the present paper x stands for the momentum fraction of the light quark 
q in the kaon (K). Taking the above two groups of DA moments for p , the parameters of 
ty p (x, k±) can be determined as, 



Group 1 : 
Group 2 : 



A\ = 382. GeV" 1 , 



A p = 422. GeV 



-i 



B\ = 0.311, C\ = 1.61, 
B 2 p = 0.257, Cl = 1.52. 



(14) 
(15) 



The distribution amplitudes for these two group of parameters are shown in Fig. (CO), where 
(pp(x,fi ) is determined by Group 1 parameters and 0p(x,yU O ) is determined by Group 2 
parameters respectively. For comparison, we also draw the distributions derived in Ref. 24, 
28 ] in Fig.([T]), i.e. ^(x, /jq) stands for the DA obtained in Ref.j^ and 4> p n (x, /iq) stands 



for that of Ref. 



]. One may observe that different from 4>p r (x, fi ) and 4> p n (x, /i ), both 
<fip(x, /i ) and (pp(x, /i ) are double humped curves and are highly suppressed in the endpoint 
region. Such feature is necessary to suppress the endpoint singularity coming from the hard- 
scattering kernel and then to derive a more reasonable results for the twist-3 contributions 
to the B — > K form factor. 

It is more convenient to transform the kaon wave functions into the compact parameter 
&K--space, which can be done with the help of the Fourier transformation 



^(x,b K ] 



d 2 k_|_ exp (— zkj_ • b^) ^(x, kj 



'[k|<l/6 K 

where \1/ stands for ^k, ^ p and ^/ a , respectively. The upper edge of the integration |kjj < 
1/bx is necessary to ensure that the wave function is soft enough 39(. After doing the Fourier 
transformation, we obtain the kaonic wave functions in the compact parameter 6^-space: 



2ttA 



K 



l + J B A -C 1 3/2 (2x-l)]exp 



x(l — x) 



-P 2 K 



x 



x 
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o — x ) 



1 — X X 



exp . , 
i) V^l-L — x ) 



ty a (x } b K ) = 1ixA a exp 

and 

= ^ [1 + B p Cl /2 (2x - 1) + C p C 2 1/ "(2a ; - 1)] exp 



^ k2± ' J (b K k ± )k ± dkj 



4/2, 



\ 1 — x a; , 



rl/6jc / £2 \ 

x / exp ( ^ 3 ] J (b K k±)k±dk± . 



o 



x(l — x) 



III. NUMERICAL CALCULATIONS 



In the numerical calculations, we adopt 

A ( ^ =4) = 250MeV, f B = 190MeV, M B = 5.279GeV, f K = 160MeV, M K = 494MeV. 

As for the phenominological parameter fijc = M^/(m s + m u ), which is a scale characterized 
by the chiral perturbation theory, we take its value to be Hk ~ 1.70 GeV. 

In the following, we first discuss the properties of F^ K (q 2 ) and F^~* K (q 2 ) that are 
calculated up to 0(l/m 2 ), i.e. to show how F+^ K (q 2 ) and F$~* K {q 2 ) are affected by the 
B-meson wave function and the kaonic wave functions. The B-meson wave functions 
and b up to next-to-leading order Fock state expansion depend on two phenomenological 
parameters A and 5. An estimate of A using QCD sum rule approach gives A = 0.57 ± 
0.07GeV [4^]. By comparing the PQCD results of the B — > tx form factor with the QCD 
LCSR results and the lattice QCD calculations, Ref. [l| shows that A = 0.55 ± 0.05GeV. As 
for the value of 5, it has been pointed out that if the contribution from the B-meson three- 
particle wave function is limited to be within ±20% of that of the WW-like wave function 
within the energy region of Q 2 G [0, ~ lOGeV 2 ], then the value of 6 should be restricted 
within the region of [0.25, 0.30] [2J. For clarity, we take the same regions as obtained from the 
B -> 7T case |l|, U for both A and 5, i.e. A G [0.50, 0.60]GeV and 5 G [0.25,0.30], to study 
the form factors F^~' K (q 2 ) and F^ K (q 2 ) in the large and intermediate energy regions. 
Furthermore, according to the discussion in the last section, the remaining uncertainty of 
the kaonic twist- 2 wave function y k is caused by the value of af(lGeV), cf. Eq.([9]). There 
we take af (lGeV) = 0.05 ± 0.02 [21f] to do our discussion. As for the twist-3 wave function 
ty p , we take two groups of parameters as shown in Eqs. fll4|ll5l) to do the calculation. 
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Next, we compare the 0(1/ m 2 ) result of the form factor with the leading order one that 
is of order 0(1/ mi) and is calculated by using the WW-like B-meson wave function, and 
also make a comparison 

energy regions. Through comparison, preferable values for the undetermined parameters can 
be found. The B — > K transition form factors F?^ K (q 2 ) and F^ K (q 2 ) have been studied 



within the framework of QCD LCSR lOj, especially at q 2 = 0, it shows 



F^ K (0) = 0.331 ± 0.041 + 0.25[af (lGeV) - 0.17], (16) 

e.g. when af (lGeV) = 0.05, F^ K (0) = 0.301 ± 0.041. More generally, F^ K (q 2 ) and 
F^^ K (q 2 ) can be parameterized in the following form [lol |: 

F^ K (q 2 ) = f as (q 2 ) + af (/i )r"(g 2 ) + af (^)f<(q 2 ) + af (»o)f a * (q 2 ), (17) 

where f as contains the contributions to the form factor from the asymptotic DA and all 
higher- twist effects from three-particle quark-quark-gluon matrix elements, f a i > a 2' a £ con- 
tains the contribution from the higher Gegenbauer term of DA that is proportional to af , 
af and af respectively. Here the factorization scale /io should be taken as 2.2GeV, since 
the functions f as > a i> a 2> a A are determined with /j, = 2.2GeV llOll . The explicit expressions 
of f as > a i> a 2> a A can be found in Table V and Table IX of Ref.jl0|. For the Gegenbauer mo- 
ments af (2.2GeV) and af (2.2GeV), we take their preferred values: af (2.2GeV) = 0.080 



and af (2.2GeV) = -0.0089 [lOj. While for af (2.2GeV), it equals to 0.793af (lGeV) with 
the help of QCD evolution. 



A. Basic properties of the form factor up to 0(l/m 2 ) 

First, we discuss the properties of F^ K (q 2 ) and F ( ^~' K (q 2 ) caused by the B-meson wave 
function. For such purpose, we fix the kaonic wave functions by setting af (lGeV) = 0.05 
and by using the Group 1 parameters for We show the B — > K transition form factors 
F^ K (q 2 ) and F^^ K (q 2 ) with 5 = 5 C = 0.275 in Fig. (T5]), where A varies within the region of 
[0.5GeV,0.6GeV]. For comparison, we show the QCD LCSR result with af (lGeV) = 0.05 
and its theoretical error (~ ±10%) [10| by a fuscous shaded band in Fig.(T5]). The results show 
that the B —> K transition form factors will decrease with the increment of A. And the best 
fit of the QCD LCSR result at q 2 = shows that A = A c = 0.525GeV. Moreover, we show 
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q 2 (GeV 2 ) q 2 (GeV 2 ) 

FIG. 2: PQCD results for the B -> K transition form factors F^ K {q 2 ) (Left) and F^~* K {q 2 ) 
(Right) with 5 = 0.275 and af r (lGeV) = 0.05. The dash-dot line, the dashed line and the dotted 
line stand for A = 0.50 GeV, 0.55 GeV and 0.60 GeV respectively. For comparison, the solid line 
comes from the QCD LCSR result as shown in Eq ■ (|17j) and the fuscous shaded band shows its 
theoretical error ±10%. 




q 2 (GeV 2 ) q 2 (GeV 2 ) 

FIG. 3: PQCD results for the B -> K transition form factors F^ K (q 2 ) (Left) and F^ K (q 2 ) 
(Right) with A = 0.525 GeV and af (lGeV) = 0.05. The dotted line, the dashed line and the 
dash-dot line stand for 5 = 0.25, 0.275 and 0.30 respectively. For comparison, the solid line comes 
from the QCD LCSR as shown in Eq. (|17|) and the fuscous shaded band shows its theoretical error 
±10%. 
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q 2 (GeV 2 ) q 2 (GeV 2 ) 



FIG. 4: PQCD results for the B -> K transition form factors F$~* K ' (q 2 ) (Left) and F^ K (q 2 ) 
(Right) with A = 0.525 GeV and 5 = 0.275. The dotted line, the dashed line and the dash-dot line 
stand for of (lGeV) = 0.03, 0.05 and 0.07 respectively. 

F+^ K {q 2 ) and F^ K (q 2 ) with A = A c = 0.525 GeV in Fig.©, where 5 varies within the 
region of [0.25, 0.30]. The results show that the B — > K transition form factors will increase 
with the increment of 5. It can be found that when setting af(lGeV) = 0.05, and by varying 
5 within the region of [0.25,0.30] and A within the region of [0.5GeV, 0.6GeV], F^ K (0) 
runs within the region of [0.23,0.34]. since the best agreement between the PQCD result 
and the QCD LCSR result at q 2 = is obtained around A c = 0.525GeV and 5 C = 0.275, we 
shall always take A = A c and 5 = 5 C to do our following calculations if not specially stated. 

Second, we discuss the properties of F^~' K (q 2 ) and F^ K (q 2 ) caused by the twist- 2 
wave function ^k, i-e. by the value of af(lGeV). For such purpose, we fix the B-meson 
wave functions by setting 5 = 5 C and A = A c and by using the Group 1 parameters for 
We show the B — > K transition form factors F^ K (q 2 ) and F^^ K (q 2 ) in Fig.(T4]) with 
af(lGeV) = 0.03, 0.05 and 0.07 respectively. It can be found that the form factors shall 
be increased with the increment of af-(lGeV), which agree with the observation of Ref . flo| . 
Furthermore, since the contribution from is sizable to that of , it is necessary to make a 
discussion on its uncertainty to the B — >• K transition form factor. Fig.© shows F^ K (q 2 ) 
and F^ K (q 2 ) with two groups of parameters for ^> p . The results are very close to each 
other due to the close shape of their <p p as shown in Fig.©, e.g. around the region of q 2 ~ 
the difference between them is less than 6%. So by taking proper transverse momentum 
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FIG. 5: PQCD results for the B -> K transition form factors F^ K (q 2 ) and F^ K {q 2 ) with 
A = 0.525GeV, 5 = 0.275 and of (lGeV) = 0.05. The dash-dot and the dashed lines are for tf p 
with Group 1 parameters Eq. (|14|) . Group 2 parameters Eq. (|15D respectively. For comparison, the 
solid lines come from the QCD LCSR with of (lGeV) = 0.05 [id]. 
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4 6 

q 2 (GeV 2 ) 



(a) (b) 
FIG. 6: PQCD results for the B -> A" transition form factor F^ K {q 2 ) with fixed A = 0.55GeF, 
5 = 0.275 and of (lGeV) = 0.05. The left dia gram is for the different kaon twist structures, ^k, 
and ^o-. The right diagram is for the different B meson structures, ^b, and A. 

dependence for the wave function where we have taken the BHL prescription for its 
transverse momentum dependence, the uncertainties from its distribution amplitude <j) p can 
be reduced. 

Finally, in order to get a deep understanding of the B — > K transition form factor, we 
discuss the contributions from different parts of the B-meson wave function or the kaon 
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wave function, correspondingly. Here we take F^ K (q 2 ) to do our discussions and the 
case of F^ K (q 2 ) can be done in a similar way. For convenience, we set A = A c , 5 = 8 C , 
af (lGeV) = 0.05 and by using the Group 1 parameters for ty p . When discussing the 
contribution from one of the kaon wave function structures, the contribution from all the 
B-meson wave function structures are summed up, and vice versa. Fig.(jS^) shows the 
contributions from the different twist structures of the kaon wave function, i.e. ^ k, 
and ^ a (the contributions from the terms involving *ff' a are included in respectively. 
One may observe that the contribution from is comparable to that of e -g- its 
contribution is about 70% of that of at q 2 ~ 0, and the contribution from ^ a is small. 
Fig.fJSb) presents the contributions from ^b-, ^ b and A respectively. It can be found that 
the contribution from ^b is about 50% — 67% in comparison to that of in the region 
of q 2 G [0, lOGeV 2 ], while the contribution from A is negligible in comparison to that of 
^b and ^/b- So the contribution from should be included for a consistent estimation 
to the next leading order. As a comparison, it can be found that under the leading order 
estimation the contribution from \I/b is only about 20% in comparison to that of \I/b at 
q 2 = 0, which agree with the rough order estimation that the contribution from is 
of order 0{\jm^. So as to the leading order estimation 0(l/mi,), the contribution from 
is usually neglected in the literature. Such difference of \I/b's contribution between the 
leading order estimation and the next-to-leading order estimation is mainly due to the fact 
that the transverse momentum dependence of the B-meson wave functions are merely a 
delta function under the WW-approximation (the leading-order estimation), while it shall 
be broadened to a certain degree according to the value of 5 by taking into account the 
3-particle Fock states' contributions (the next-to- leading order estimation), cf. fig. (2) of 
Ref.Q]. So qualitatively, the contributions from ^b shall be raised to a certain degree for 
the next-to-leading order case, due to the less suppression of the end-point region (£ — > 0) 
from the transverse momentum distributions than that of the leading order case. And 
then the naive order estimation for the contribution of \I/b is no longer correct, and the 
contributions from \1/ b and ^ b are both important in the next-to- leading order calculation. 
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B. Comparison with the leading order results 



The WW-like B- meson wave functions in the compact parameter fen-space can be found 
in Ref.[l|. Taking the WW-like wave functions and cutting off the power suppressed terms 
in the hard scattering amplitude, we can obtain the leading order results (d?(l/m&)) for the 
form factors F^ K (q 2 ) and F B "* K {q 2 ). Strictly, one should cut off the contribution from \1/ r 
to obtain the leading order estimation, since is power suppressed in comparison to ^b- 
However for easy comparison with the results in the literature, e.g. Ref . [7| , we keep \I/ b in the 
leading order estimation. For convenience, we take A = A c , 5 = 5 C , af(lGeV) = 0.05 and 
by using the Group 1 parameters for the wave function ^ p to do a comparison of the leading 
order results with the total results that include the contributions up to order 0(l/m 2 ). It 
can be found that the leading order results are smaller than the total results by about 25% in 
the large recoil region, e.g. at q 2 = 0, the leading order F^ K (0) = 0.229. One may observe 
that a larger leading order estimation has been obtained in Ref. [7|, which shows F+~^ K (0) = 
0.321 ± 0.036. We argue that the present leading order estimation is more reliable, and the 
larger value of F+~^ K (Q) derived in Ref.|7| is mainly due to the following two reasons: 1) Even 
though the Sudakov and threshold resummation factors shall kill the endpoint singularity 
of the process 5], ?| Q, IlS], the transverse momentum dependence of kaonic wave functions 
are still important to give a more reliable PQCD estimation, which is similar to the cases 

ft rn 

of B — > 7r form factor |l| and the pion electromagnetic form factor [29j. In Ref. [7j] the 

transverse momentum dependence of kaonic wave functions are lacking, i.e. the distribution 

amplitude other than the wave function is used. While in our present calculation, the BLTL- 

prescription is adopted for the kaonic transverse momentum dependence. As for the wave 

function ^f p (x, kjj, such transverse momentum dependence will results in a double humped 

DA p as shown in Fig. ([1]) and then it shall give more effective suppression in the end-point 

Q 

region than the one used in Ref.[7j. In fact, it can be found that the contributions from the 
end point region shall always be overestimated without taking the transverse momentum into 
the twist-3 wave function ^ p (x, kj_) 3 . Furthermore, by taking proper transverse momentum 
dependence for the wave function the uncertainties from its distribution amplitude 4> p 
can be reduced as has been discussed in Sec. III. A; 2) the distribution amplitude of with 



For example, a detailed discussion on the model dependence of pionic twist-3 wave function ^ p (x,]ij_) 
can be found in Ref. 291. 
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a much bigger value of af^lGeV), i.e. af-(lGeV) = 0.17, is adopted by Ref.j?]. Since the 
form factors increases with the increment of af (lGeV), a larger value of af (lGeV) shall 
increase the form factors. 

Furthermore, by varying A within the region of [0.50, 0.55], the uncertainty caused by A 
is the biggest and is of order (1/rrib). While by varying 5 within the region of [0.25,0.30], 
the uncertainty caused by 5 is smaller and are of order (1/ml). This can be qualitatively 
explained as that A is the characteristic parameter that determines the leading Fock-state 
behavior of the .B-meson wave functions, while 5 is the characteristic parameter that deter- 
mines the higher Fock-state's behavior of the .B-meson wave functions. The uncertainties 
from af and ^ k are less than 10% in the large recoil region. 



C. Comparison with the LCSR results 



The B — > K transition form factor have been analyzed by several groups under the QCD 



LCSR approach 



as 
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111 ]. New sum rule for B — > K is derived from the correlation 



in which the 



functions expanded near the light cone in terms of B-meson distributions 
contributions of the quark-antiquark and quark-antiquark-gluon components in the B-meson 
are taken into account. It has been found that the B — > K transition form factor in the 
large recoil region does not receive contributions from the 3-particle B-meson DA's. One may 
observe that if substituting the B-meson DAs, which are derived by doing the integration 



over bs in Eqs. (|4"|5]) . into the formulae of Ref. 



then one can obtain the same results as that 



of Ref. ja , since our B-meson DAs are close to the exponential model wave functions adopted 
in Ref.[8|. Furthermore, one may observe that the result of F^ K (0) = 0.31 ± 0.04 under 
the condition of af (lGeV) = 0.05 ± 0.03 agrees well with our present PQCD estimation. 
Secondly, a systematic QCD LCSR calculation has been done in Ref.floT] by including the 
one-loop radiative corrections to the twist- 2 and twist-3 contributions, and leading-order 
twist-4 corrections. Some comparison of their results with our present one can be found in 
Figs. (l2l[3l4l) . which also shows a good agreement within reasonable errors. For example, from 



Eq. (1171) it can be found that the uncertainty of form factor caused by af (lGeV) within the 
region of [0.03, 0.07] is less than 5%, which is consistent with our present result as shown in 
Fig.©. 



18 



IV. DISCUSSION AND SUMMARY 



In this paper, we have examined the B — > K transition form factor in the PQCD approach 
up to order 0(1/ ml), where the transverse momentum dependence for the wave function, the 
Sudakov effects and the threshold effects are included to regulate the endpoint singularity 
and to derive a more reasonable result. We have confirmed that the PQCD approach can 
be applied to calculate the B — ► K transition form factor in the large recoil regions. We 
emphasize that the transverse momentum dependence for both the B meson and the kaon 
is important to give a better understanding of the B —> K transition form factor. Fig.([6K) 
shows that the contribution from the pionic twist-3 wave function ty p is sizable in comparison 
to that of ^ Ki and the contribution from is small. While Fig.([Hb) shows that by using the 
B-meson wave functions up to next-to-leading order Fock state expansion, the contributions 
from and 1 5r are important. 

In Refs.jl, J, we have shown that the results from the PQCD approach, the lattice 
QCD approach and the QCD LCSRs are complementary to each other and by combining 
the results of those three approaches, one can obtain an understanding of the B —>■ ix 
transition form factor in the whole physical regions. And the best fit of the PQCD results 
with that of the QCD LCSR results in the large recoil region can be obtained by taking 
A G [0.50,0.60] and 5 G [0.25,0.30] 2j. In the present paper, we show that within the 
regions of A G [0.50,0.55], S G [0.25,0.30] and af(lGeV) G [0.03,0.07], the PQCD results 
on the B — > K form factor in the large recoil region also agree well with that of the QCD 



LCSR results 



101 ] . Our present PQCD results in some sense is more reliable than the 



LCSR calculations due to the fact that by taking the transverse momentum dependence 
properly for the wave functions the soft endpoint singularity have been effectively suppressed, 
e.g. as is shown in Sec. III. A the difference between the two models for ty p is less than 
6% in the large recoil region, while for the LCSR approach large uncertainty comes from 
the kaonic twist-3 DA <p p that is not too well-known. By running the parameters within 
the above regions, we obtain F^ K {0) = 0.30 ± 0.04. Finally, to illustrate the SU f {3)- 
breaking effects, we calculated the ratio with the help of the B — > tx results in Ref.j^]: 
[F+fl K {ty/F+jo*( Q )] = 1-13 ± 0.02, which favors a small SU /(3)-breaking effects. 
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